Observation of the Kohn anomaly near the K point of bilayer graphene 
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The dispersion of electrons and phonons near the K point of bilayer graphene was investigated 
in a resonant Raman study using different laser excitation energies in the near infrared and visible 
range. The electronic structure was analyzed within the tight-binding approximation, and the 
Slonczewski-Weiss-McClure (SWM) parameters were obtained from the analysis of the dispersive 
behavior of the Raman features. A softening of the phonon branches was observed near the K 
point, and results evidence the Kohn anomaly and the importance of considering electron-phonon 
and electron-electron interactions to correctly describe the phonon dispersion in graphene systems. 
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Graphene systems exhibit a strong electron-phonon 
coupling at special points in the Brillouin zone, that 
softens the phonon energy and gives rise to kinks in 
the phonon dispersion (infinities in duj/dq), which are 
called Kohn anomaly [l], Q ■ This effect has been demon- 
strated experimentally at the T point of monolayer and 
bilayer graphene using gated Raman scattering exper- 
iments [3, 0, [H, H 0, HI • However, the electron-phonon 
coupling is expected to be stronger at the K point [ij , but 
Raman experiments in graphene systems performed with 
visible light cannot probe phonons near the K point 
Moreover, in the case of AB-stacked bilayer graphene, 
due to its special electronic and phonon structure [9j , Ra- 
man experiments involve phonons closer to the K point 
when compared to monolayer graphene [Io| . This work 
presents a resonance Raman investigation of AB-stacked 
bilayer graphene using many laser lines in the near- 
infrared (near-IR) and visible range. The Kohn anomaly 
for both symmetric (S) and anti-symmetric (AS) phonons 
was evidenced, and results show the importance of con- 
sidering higher renormalization terms such as electron- 
electron interactions to correctly describe the phonon dis- 
persion near the K point [ll[ . These effects are especially 
relevant for understanding the dispersion of electrons and 
phonons [lp. [l3l [Llj and transport properties in this novel 
material [9j. 

In a previous resonance Raman study of bilayer 
graphene performed in the visible range [la ], the elec- 
tronic structure of bilayer graphene was probed by ana- 
lyzing the dispersion of G' Raman band (also called 2D 
band) as a function of the laser energy, and described 
within a tight-binding approximation [lq . UtI [l8| by de- 
termining the nearest-neighbors hopping parameters 70, 
71, 73 and 74. It was shown in this work that a linear 
iTO phonon dispersion provided a good fit of the exper- 
imental data obtained with visible photons [l5| . 

In the present work, we have extended the range of 
laser energies, measuring the G' Raman band with many 
laser lines in the range 1.33 to 2.81 eV. The measurements 



in the near-IR range (1.33 to 1.69 eV) are especially rele- 
vant since we can probe phonons that are much closer to 
the K point. The analysis of the low energy data allowed 
us to observe a non-linear softening of the phonon branch 
near the K point, and the significant splitting of the sym- 
metric (S) and anti-symmetric (AS) phonon branches. In 
particular, we show that the phonon softening is stronger 
for the S branch. Concerning the electronic structure, we 
have also considered the in-plane second-neighbor hop- 
ping parameter, which is expected to be of the same or- 
der as the out-of-plane nearest-neighbor parameters, to 
describe the G' Raman band dispersion. 

The graphene samples used in this experiment were 
obtained by a micro-mechanical exfoliation of graphite 
(Nacional de Graffte, Brazil) on the surface of a Si sam- 
ple with a 300 nm Si02 coverage. The laser power was 
kept at 1 mW in order to avoid sample heating. We used 
He-Cd, Ar-Kr and dye lasers for the laser lines in the vis- 
ible range (1.91-2.81 eV) and a TkSaphire laser for the 
excitation in the near-IR range (1.33-1.69eV). 

Figure [1] shows the G' band of bilayer graphene 
recorded with 19 different laser lines between 1.33 and 
2.81 eV (932 to 440 nm). We can see that both the po- 
sition and the shape of the G' band are strongly depen- 
dent on the energy of the exciting laser. The G' band in 
graphene systems comes from an intervalley double reso- 
nance (DR) Raman process [l^, [2(| that involves one ini- 
tial electronic state with wavevector k near the K point, 
one intermediate electronic state with wavevector kl near 
the K' point, and two in-plane transverse optical (iTO) 
phonons with wavevectors q — k + k' [lOj] . Since photons 
with different energies excite electrons and phonons with 
different wavevectors k and q, respectively, the dispersion 
of electrons and phonons near the K point can be mea- 
sured in when the energy of the incident photons can be 
tuned [l(j. 

In the case of AB-stacked bilayer graphene, the elec- 
tronic branches split into two valence bands [tv\ and 7^) 
and two conduction bands (n* and tt|) The iTO 
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FIG. 1: The Raman G' band of bilayer graphene recorded 
with 19 different laser lines between 1.33 and 2.81 eV (932 to 
440 nm). 

phonon branch also splits into two branches, related to 
the symmetric (S) and anti-symmetric (AS) phonons. 
It has been shown that the DR scattering process oc- 
curs preferentially in the high symmetric TKM direc- 
tion [U [22j , involving electrons in the KT direction and 
phonons in the KM direction. 

Group theory analysis for bilayer graphene [23| pre- 
dicts four distinct DR processes (Pn, P22, P12 and 
P21) along the TKM direction, which are illustrated in 
Figs. H)^ a) and (b). Since the lower and upper conduction 
bands (n^ and 7rJ) belong to different irreducible repre- 
sentations, the S phonons (Tj symmetry) are associated 
with the Pn and P22 processes [see Fig.^a)] involving 
electrons with same symmetry {tt\ — > 7r* or 7^ — > tt|) 
whereas the AS phonons (T2 symmetry) occur for pro- 
cesses P12 and P21 [Fig.[2[b)] involving electrons with 
different symmetries (n* — > 7r| or tt^ — > 7rJ ) [l]| . 

Each of these processes (Pn, P22, P12 an d P21) is 
responsible for one peak in the G' band of AB-stacked 
bilayer graphene. Fig.[2]Jc) shows the Raman spectrum 
of the G' band measured with 1.57 eV laser energy and 
fitted with four peaks. All have the same FWHM of 
24 cm -1 , which is the FWHM of the single band in mono- 
layer graphene [TH [Hj]. We see in Fig^a) that the 
Pn process involves the phonon with higher wavevec- 
tor, while the phonon with smaller wavevector gives rise 
to the P22 process. Therefore, the P22 process is espe- 
cially relevant since it involves phonons closer to the K 
point, even when compared to the phonon in monolayer 
graphene probed by the same laser line. Since the iTO 
phonon frequency increases with increasing q, we con- 
clude that the lowest and highest frequency peaks of the 
G' band are associated with the Pn and P22 processes, 
respectively. The two intermediate Raman peaks are due 
to the P12 and P21 processes, as shown in Figj2jc). 

All G' bands shown in FigQ] were fitted by four 
Lorentzian curves, and laser energy dependence of the 
Lorentzian peak positions are plotted in Fig(3Ja-d). In 
order to analyze the experimental dispersion of the G' 




FIG. 2: (Color online) (a) Pn and P22 DR Raman processes 
involving the symmetric phonon, (b) P12 and P21 DR Ra- 
man processes involving the anti-symmetric phonon (c) G' 
Raman band of bilayer graphene measured with 1.57 eV laser 
energy and fitted with four Lorentzian peaks with FWHM of 
24 cm -1 , (d) Bilayer graphene atomic structure showing the 
hopping processes associated with first-neighbor interactions 
(70, 7i, 73 and 74), and second neighbor interaction t' . 



peaks shown in Fig[3]Ja-d), we need to consider the dis- 
persion of both electrons and phonons near the K point. 
The electronic dispersion will be analyzed here using the 
tight-binding approximation, which was first introduced 
by Wallace 1181 . using the Slonczewski-Weiss-McClure 
(SWM) [l6|, Il7l | model for graphite, and the phonon dis- 
persion is obtained from the fit of the experimental data. 

In the previous resonance Raman study in bilayer 
graphene performed in the visible range [la ], a linear 
phonon dispersion was considered to fit the G' peak posi- 
tions versus E; aser data. Fig. [3Ja) shows that the fitting 
of the data of the present study considering the linear 
phonon dispersion and the SWM parameters used in ref- 
erence 15] fails for the experimental points in the near-IR 
region and, in particular, for the data associated with the 
P22 process involving phonons closer to the K point. 

In order to fit the low energy experimental data in Fig. 
[31 we have considered a non-linear relation for the iTO 
phonon dispersion, given by a second-order polynomial 
w(q) = A + Bq + Cq 2 . FiglJJb) shows the fit consid- 
ering the same non-linear phonon dispersion for the S 
and AS phonon branches, and the 70, 71, 73, 74 parame- 
ters. FigJJJc) shows the fit using the same SWM parame- 
ters as in Fig[3jb) but considering two distinct non-linear 
phonon dispersions for the S and AS branches. As we 
can see in Fig[3Jb) and (c), the fitting in the low energy 
range is improved considering the non-linear dispersion, 
but different dispersions for the S and AS branches are 
needed to obtain a good fit of the experimental data. All 
fitting parameters are shown in Table [TJ 

In the fittings shown in Figs. [3fa)-(c), we have consid- 
ered only the first-neighbor parameters 70, 71, 73, 74. In 
principle, we could also introduce higher-order terms and, 
in particular, the in-plane second-neighbor parameter t' , 
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FIG. 3: Laser energy dependence of the peaks of the G' bands, 
fitted considering four different approximations, (a) Fitting 
considering linear phonon dispersion and the SWM parame- 
ters 70, 71, 73, 74 (same approximation as in Ref. [lap, (b) 
Fitting considering the same non-linear phonon dispersion for 
the S and AS branches, and the SWM parameters 70, 71, 73, 
74. (c) Fitting considering different non-linear phonon dis- 
persions for the S and AS phonon branches, and the same 
SWM as in (b). (d) Fitting using te same phonon dispersions 
as in (c), the first-neighbor SWM parameters 70, 71, 73, 74, 
second-neighbor in-plane t' parameter, and the A parameter, 
which represents the difference in energy of the sublattices A 
and B. All fitting parameters used in approximations (a), (b), 
(c) and (d) are shown in Table [I] 



which is expected to be of the same order of magnitude as 
the out-of-plane first-neighbor parameters. Figure [3^d) 
shows the fit of the experimental data considering also t' 
and A, which represents the difference in energy of the 
sublattices A and B. As we observe in Figure G^d), the 
resultant fit is slightly improved, mainly due to the use 
of larger number of fitting parameters. 

It is important to emphasize that slightly different val- 
ues of the 7 parameters are found when we include t' and 
A. A good fit can always be obtained for 70 values rang- 
ing between 2.9 and 3.1 eV. Concerning the 71 parameter, 
reasonable fits could only be obtained for 71 < 0.35 eV, in 
disagreement with values of 71 up to 0.4 eV proposed in 
the literature [H, [2(| [13] • Considering the 73 parameter, 
the best fit is obtained when 73 w 0.1 eV, and a reason- 
able fit cannot be obtained for values of 73 > 0.15 eV. 
Once again, this value is smaller than others found in the 
graphite literature (73 « 0.30 eV [H, H| [H). Recent 
infrared studies in exfoliated bilayer graphene consider 
73 w 0.30 eV [U HH, but this value is not extracted 
directly from the experiments [25|, [27|. Notice that 73 
is related to the trigonal warping effect (TWE) at very 
low energies and gives rise to electron-hole pockets [33J 
in the energy scale of ~ 2 meV, which is not accessible 
in Raman experiments. However as the energy increases, 
there is also a TWE due to the symmetry of the lattice, 
which is probed by our experiment. These two effects 



can explain the different 73 values found with distinct 
experimental techniques. 

If we consider only nearest neighbor parameters to de- 
scribe the electronic structure, the best fit is obtained for 
74 « 0.15 eV, which is in close agreement with our pre- 
vious experiment [HI and with Refs. [25|, [2(|. However, 
smaller values of 74 provide a good fit when the second- 
neighbor t' parameter is included. In fact, both 74 and t' 
parameters are associated with the asymmetry between 
electrons and holes in bilayer graphene. Finally, reason- 
able fits can be obtained for different small positive and 
negative values of A (|A| < O.OleV). 

The analysis of the experimental data also gives us 
the phonon dispersion for the symmetric (S) and anti- 
symmetric (AS) iTO phonon branches near the K point. 
As shown in Figs.[3^c) and (d), two distinct phonon 
branches, associated with symmetric and anti-symmetric 
iTO phonons, are necessary to fit the experimental data. 
The values of the polynomial parameters A, B and C for 
the S and AS branches are shown in Table|TJ We observe 
that the quadratic coefficients C assume negative values, 
showing that the slope of the phonon dispersion increases 
with decreasing q values. This is direct evidence of the 
Kohn anomaly for the iTO phonon branches, which are 
expected to exhibit a kink at the K point. 

The experimental data shown in Fig. [3] can be di- 
rectly plotted in a phonon energy dispersion relation 
(cOph versus q) by eliminating the laser energy (Ei aser ) 
in the double resonance conditions [l5| . Fig. [J] shows 
the phonon dispersion of the S (full circles) and AS 
(open circles) branches obtained from the resonance Ra- 
man results, considering the same TB parameters as in 
FiglUc). Fig. |4] also shows the theoretical phonon dis- 
persion near the Dirac point of the iTO phonon branch 
of monolayer graphene calculated using the tight-binding 
approximation by Popov et a l. 1 341] (solid curve) and us- 
ing DFT by Lazzeri et al. [Tl| (dashed curve) within 
the GW approximation, where electron-electron interac- 
tion is taken into account. Notice that the phonon dis- 
persion calculated by tight-binding fails to describe the 
data for lower phonon energies, which are in good agree- 
ment with the calculations within the GW approxima- 
tion. This result shows the importance of considering 
electron-phonon and electron-electron interactions in or- 
der to correctly describe the phonon dispersion near the 
K point of graphene systems. 

Another interesting observation is the fact that the 
phonon softening is stronger for the S phonon branch 
as shown in Fig. [4] This result is in agreement with 
the calculation performed by Ando et al. [35| . which 
predicts a stronger phonon renormalization for the zone- 
center symmetric phonon, due to distinct selection rules 
for interaction of S and AS phonons with intra- valley and 
inter- valley electron-hole pairs. The same type of selec- 
tion rule is expected to occur for phonons near the K 
point. 

In summary, the dispersion of electrons and phonons 
of bilayer graphene was investigated by performing a res- 
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TABLE I: Values of the SWM parameters (in units of eV) and the iTO phonon dispersion parameters (w(q) = A + Bq + Cq 2 ) 
obtained from the four different fits of the experimental data shown in Fig[3](a), (b), (c) and (d). 
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FIG. 4: (Color online) Experimental iTO phonon dispersion 
of bilayer graphene near the K point of the S (full circles) 
and AS (open circles) branches. The solid and dashed curves 
correspond, respectively, to the theoretical iTO phonon dis- 
persion of monolayer graphene calculated using tight-binding 
by Popov et al. [34| and using DFT by Lazzeri et al. [ll| 
within the GW approximation. 



onance Raman study of the the G' Raman band, using 



laser energies from the visible to the near-IR range. The 
measurements in the near-IR range are especially rele- 
vant since we can probe phonons that are much closer 
to the K point. The electronic structure was analyzed 
within the tight-binding approximation, considering first 
and second neighbor interactions. We have obtained ac- 
curate experimental data for the phonon dispersion of the 
iTO branches near the K point, and the phonon branch 
softening reveals the K point Kohn anomaly in bilayer 
graphene. We have shown that the phonon renormal- 
ization is stronger for the S phonon branch. Our results 
agree with the phonon dispersion calculation which takes 
into account electron-electron interaction in graphene 
systems, which plays an important role to correctly de- 
scribe the Kohn anomaly near the K point. 
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